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Abstract 


We review one of the most fruitful areas in cosmology today that bridge 
theory and data — the temporal growth of large-scale structure. We go 
over the growth’s physical foundations, and derive its behavior in sim- 
ple cosmological models. While doing so, we explain how measurements 
of growth can be used to understand theory. We then review how some 
of the most mature cosmological probes — galaxy clustering, gravita- 
tional lensing, the abundance of clusters of galaxies, cosmic velocities, 
and cosmic microwave background — can be used to probe the growth 
of structure. We report the current constraints on growth, which are 
summarized as measurements of the parameter combination fog as a 
function of redshift, or else as the mass fluctuation amplitude parame- 
ter Sg. We finally illustrate several statistical approaches, ranging from 
the “growth index” parameterization to more general comparisons of 
growth and geometry, that can sharply test the standard cosmological 
model and indicate the presence of modifications to general relativity. 
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1 Introduction 


In the standard cosmological model, the seeds of structure laid out by inflation 
were amplified by the influence of gravity. Initial perturbations of the size 
ô ~ 1075, where 


ô 

= (1) 
p 

is the density perturbation in matter, has been amplified by gravity over the 


billions of years that followed the inflationary epoch. The general fact that 


ô 
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structure grows over time is a very well-established feature of the hot big- 
bang cosmological model. The growth of structure can be measured through 
observations of positions and motions of galaxies as a function of cosmic time, 
and by comparing the inferred overdensities to the initial conditions, ô ~ 1075, 
that can be observed in the temperature anisotropies of the cosmic microwave 
background (CMB). 

A cosmological model predicts the temporal growth of structure. Such a 
prediction is a key ingredient in theoretical calculations of the various observ- 
able quantities that can be experimentally measured. For example, weak 
lensing shear integrates the power spectrum of structure along the line of sight 
(as we discuss in Sec. 3.2), which is in turn is determined by the growth of 
structure. As another example, the amplitude of galaxies’ peculiar velocities, 
as well as the spatial correlations of galaxy positions, are both largely deter- 
mined by the growth of structure. In both of these examples, a larger rate 
of growth corresponds to a larger signal (assuming fixed initial conditions). 
Turning the argument around, measurements of the large-scale structure and 
inference of its statistical properties inform us about the growth of structure. 
In turn, constraints on the growth of structure help pin down the parameters 
of the cosmological model, and are especially sensitive to the properties of dark 
energy. 

Over the past two decades, the advent of massive new datasets from 
large-scale structure has brought into focus the importance of the growth 
of structure. This is principally because growth is a sensitive probe of both 
dark energy and modified gravity, as hence its measurements can distinguish 
between the two. Specifically, in the context of general relativity, growth is 
mathematically related to the geometrical measures (such as distances), and 
knowing the former determines the latter or vice versa. However, this relation 
between growth and geometry that assumes standard relativity is generally 
broken when gravity is modified. Comparing growth to geometry thus enables 
stringent tests of modify gravity. 

The purpose of this review is to present a rather high-level review of the 
growth of structure, aimed at a non-expert astrophysicist. Our main goal to 
lay out the theoretical foundations for basic studies of growth. We also wish 
to illustrate, at a basic yet quantitative level, how growth enters the various 
cosmological probes, and how they can be used to measure growth. Finally, 
we review some of the most promising strategies to use growth to better 
understand dark energy and modified gravity. This review focuses on the big- 
picture foundations of the subject and the cosmological utility of the growth 
of structure (as well as current constraints), and updates and complements the 
previous major review of the subject [1] which focused on future surveys and 
expected systematic errors in the associated cosmological probes. 

This review is organized as follows. In Sec. 2 we discuss the theoretical 
background useful for understanding growth of structure. in Sec. 3, we review 
the observations that are sensitive to growth, and compile recent constraints 
on the growth of structure. In Sec. 4, we discuss ways in which the growth of 
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structure can be used to probe dark energy and general relativity. We conclude 
in Sec. 5. 


2 Theoretical background 


Here we outline the theoretical predictions for the growth of cosmic structure. 
We start with density perturbations, introduce the (linear) growth of structure, 
derive an equation that governs its evolution in time, and study its solutions. 


2.1 Temporal evolution of linear density perturbations 


In our theoretical treatment, we will specialize in linear theory, and assume 
small fluctuations, with || « 1. Following a standard, non-relativistic, 
perturbation-theory analysis that combines the continuity, Euler, and Poisson 
equations, one arrives at (e.g. [2]) 


07 ők a ðk kee 2T 2 
PA + 25 = (4ncpar na ) Or 3 qa" Sk, (2) 


where T is temperature, a is the scale factor, c, is the speed of sound, and 
pm is the mean matter density; all of these quantities depend on cosmic time 
t. Further, ôk is the mode in the Fourier expansion of the overdensity field at 
some wavenumber k 


y(t) = Se | se ner, (3) 


Here V is the volume of the larger region over which the perturbations are 
assumed to be periodic (note that k and r are both comoving quantities). 
Similarly, Sk is the Fourier mode of entropy perturbations. 

We further specialize in isentropic initial conditions where there is no fluc- 
tuation in entropy in the initial conditions, so that Sk « VS = 0. [Confusingly, 
these fluctuations are most often called adiabatic initial conditions, which 
strictly speaking implies S= 0, not VS = 0; we henceforth adopt this impre- 
cise but popular nomenclature.] In the presence of adiabatic initial conditions, 
fluctuations in various components (matter, radiation, neutrinos etc) are pro- 
portional to each other, and the overall curvature fluctuation is nonzero. This 
is the kind of initial condition that inflation typically predicts, and that cur- 
rent data favor. Finally, we specialize in fluctuations on scales larger than the 
Jeans scale, sok < kyjeans = V/4tGpm (a/c). This is a reasonable assump- 
tion since, after recombination, the Jeans scale corresponds to mass of order 
10°Mo, which is much smaller than the structures that we will be interested in 
which are roughly mass of a galaxy or larger, so > 10!?Mo. In this regime, we 
can drop the term that is proportional to k?. With all that, Eq. (2) simplifies 
to ra 

6 + 2H6 — 4nGpmô = 0. (4) 
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Notice that we have dropped the subscript , in the overdensity, as there is no 
more wavenumber dependence in this equation. 

Equation (4) describes the evolution of the density perturbations on scales 
0.001 hMpc™}! < k < 0.1hMpc™!, which is where the most usable observa- 
tional data reside. This is a second-order ordinary differential equation for the 
matter overdensity 6. We will solve it below to get solutions for the linear 
growth of structure in simple cosmological models. Let us first recapitulate the 
assumptions that were assumed in deriving this equation: 


e General theory of relativity; 


Adiabatic initial conditions; 


Sub-horizon scales (i.e. relativistic effects are ignored); k >> Hp; 


Spatial scales above the Jeans length (k < ky); 


Linear theory (ô < 1), which corresponds roughly to scales k < 0.1 h Mpc7! 
today. 


The breakdown of any of these assumptions could lead to the invalidation 
of Eq. (4), which would be manifested as a difference between the predicted 
and observed growth. The growth of structure can of course be theoretically 
predicted when either one of the assumptions above is relaxed, but those pre- 
dictions tend to be more detail-dependent, have functional or numerical forms 
that are necessarily more complicated, and are also in the regimes where the 
comparison to data is more difficult (see the discussion in the bulleted list 
just below). Therefore, we will largely stick to studying the growth as given in 
Eq. (4), with the assumptions given in the bullet points above. Before proceed- 
ing, however, we wish to comment further about a few of these assumptions 
and situations in which relaxing them may be relevant: 


e Beyond general relativity: Our assumption of general relativity in 
deriving Eq. (4) deserves particular attention, as modified-gravity theories 
generically lead to the scale dependence of 6(t). Observing this scale depen- 
dence, and ruling out the systematic errors as the cause, would present 
striking evidence for the presence of modified gravity. There has been 
growing interest in trying to constrain the scale-dependent predictions of 
the growth of structure with current and future surveys [3-8]. 


e Near-horizon scale: Equation (4) is modified on scales approaching the 
Hubble distance, k ~ Hp œ 107° h Mpc™t. This is due to relativistic effects 
which can be precisely quantified [9-14]. Detecting these relativistic effects 
would be a very interesting test of the standard cosmological model, is chal- 
lenging given that they appear on very large scales, but may nevertheless 
be possible with forthcoming large-scale structure surveys [15-20]. 


e Beyond linear theory: this assumption is the most “ready” to be relaxed, 
as a large portion of current observations lies in the quasi-linear or fully 
nonlinear regime. Current surveys (e.g. KiDS [21], DES [22], and HSC [23]) 
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explicitly throw out scales that are nonlinear, but even then some infor- 
mation comes from quasi-linear regime where Eq. (4) does not hold. This 
puts premium on our ability to model the growth of structure in this regime 
with a combination of numerical and analytic tools. Such predictions are 
available in ACDM (cosmological model with vacuum energy and cold dark 
matter) [24-34], but more recent work has extended this to dark energy 
models with the equation of state parameterized by (wo, Wa) [85-39], as well 
as modified-gravity models [40-51]. 


2.2 Growth in simple cosmological models 


The growth equation (4) can of course be solved numerically for an arbi- 
trary cosmological model — the input required is the scaling of the Hubble 
parameter with time, H(t), and that of the matter density, pm (t). However, 
it is instructive to derive scalings in single-component universes, where the 
expansion is dominated by a single component — matter, radiation, or dark 
energy. 

Let us start with the radiation-dominated case. Then a œx t!/?, so that 
H(t) = a/a x 1/(2t). Also, note that the last term in Eq. (4) is negligible 
because the Hubble parameter is dominated by the radiation and not matter 
density, so that 4rGpm < H?. Therefore, we need to solve the equation 6 + 
2H = 0, with H = 1/(2t). Its solution is 


ô(t) = Ay + Ag Int (radiation dominated), (5) 


where A, and Ag are some constants. Thus, in the radiation-dominated regime, 
the fluctuations grow only very slowly — logarithmically with time. 

In the flat, matter-dominated (Einstein-de Sitter) case, a « t?/3 so that 
H(t) = a/a x 2/(3t), while the Hubble parameter is dominated by matter 
density, so that 4rGpm = (3/2)H?. Let us assume that ô(t) œ t”; then the 
growth equation simplifies to n(n — 1) + $n — 2 = 0. Its solutions are easy to 
obtain: n = +2/3 and —1. Hence 


6(t) = Byt?/? + Bot"! (matter dominated), (6) 


where Bı and By are some constants. Since a(t) œx t?/3 in the matter- 
dominated era, the growing mode of the perturbations grows proportionally to 
the scale factor, 5(t) « a(t) œ t?/3. This scaling is of the utmost importance, 
as the universe spends of order 10 billion years in the matter-dominated era 
— from the matter-radiation equality 50,000 years after the Big Bang, to the 
onset of dark energy a few billion years ago. During that time, structures in 
the universe grow appreciably, all thanks to the 6 « a scaling. 

Finally, in the dark-energy-dominated era, which will presumably take 
place in the future when dark energy dominates completely, the scale factor 
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grows exponentially', a x e#', so that H(t) = H, = const. [Here we are rep- 
resenting dark energy by a ave model, that of the cosmological constant A, 
which mathematically represents vacuum energy which has a constant energy 
density across cosmic time.] Also, note that the last term in Eq. (4) is neg- 
ligible since the matter density is negligible relative to vacuum energy in H. 
Therefore, we need to solve the equation 6 + 2H, 6 = 0 whose solution is 


S(t) = Cy + Coe? ~ ~ const (Lambda dominated), (7) 


where C and C2 are some constants, and where the exponentially decaying 
term becomes negligible quickly. Therefore, density perturbations do not grow 
at all in a Lambda-dominated universe. We are all witnesses to this effect 
today, as our universe with 70% dark-energy and 30% matter displays severely 
suppressed structure formation (relative to the matter-only scenario), some- 
thing that is readily observed in the cosmological data that probe z < 1 using 
tests that we describe in Sec. 3. 

Summarizing, the temporal evolution of linear density perturbations is 
simple, single-component cosmological models is: 


2/3 xa (matter dominated) 
d(a) x 4 In(t) « In(a) (radiation dominated) (8) 
const (Lambda dominated), 


where a is the scale factor. 


2.3 Dimensionless linear growth function 


It is useful to cast Eq. (4) in a dimensionless form. Let us introduce the linear 
growth function 
D(a) = H(A)’ (9) 


(1 
where D at the present time is unity, D(1) = 1. The original equation deter- 
mining the linear growth of structure, Eq. (4), can now be set in a dimensionless 
form as 


d?g 
dln a? 


dg 
dlna 


2 


+ [5 — 3w(a)Npe(a)] +3[1 — w(a)| Qne(a)g = 0, (10) 
where g(a) is the growth suppression factor — that is, growth relative to that 
in Einstein-de Sitter (Qm = 1) universe. The linear growth suppression factor 
g is related to the growth factor D implicitly via 


paoa O (11) 


1 The same happens during inflation, when the universe is completely vacuum-energy dominated. 
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Fig. 1 Illustration of the linear growth function D(z) and the growth suppression factor 
g(z), as a function of redshift z, in the standard flat ACDM cosmological model. For com- 
parison, we also show the results in the flat model with dark matter only (Qw = 1), the 
Einstein-de Sitter (EdS) model. Either D(z) or g(z) fully describes the linear growth of 
fluctuations; see text for details. 


Note that the linear growth function depends on the Hubble parameter H(a), 
as well as the matter density Qm. Thus the linear growth function depends 
on basic cosmological parameters; in the flat ACDM model for example, 
D = D(a,Qm), and same for g. In particular, the growth function does not 
depend on wavenumber k. [This is no longer true in modified-gravity models as 
explained above, but also in models with massive neutrinos, which introduce 
a mild scale-dependence to the linear growth function [52]. 

Figure 1 shows the temporal evolution of D(z) and g(z) in the standard 
cosmological ACDM model with Qm = 1—Q,a = 0.3 (in what follows, we 
will frequently switch between redshift z and scale factor a = 1/(1 + z)). For 
comparison, we also show the results in the flat model with dark matter only 
(Qm = 1), the Einstein-de Sitter (EdS) model. In ACDM, the onset of dark 
energy at late times (z < 1) causes these functions to deviate from their EdS 
behavior (D(a) = a or D(z) = 1/(1+ z), and g(z) = 1). The late-time decrease 
in either g(z) or D(z) fully specifies the suppression in the growth of structure 
(in linear theory) due to dark energy. In our fiducial flat ACDM cosmological 
model with Qm = 1 — Qpe = 0.3, the present-day value of the suppression 
factor is g(z = 0) ~ 0.78. 

Linear growth function enters the linear matter power spectrum in a 
straightforward way. The matter power spectrum, P(k), quantifies the amount 
of structure — the power — at each wavenumber k. The power spectrum 
is defined in terms of the two-point correlation function of Fourier-space 
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overdensity Od, as 

(ôk div) = (27)? 5 (k — k’) P(k) (12) 
where the angular brackets denote ensemble average, and 5“) is the Dirac delta 
function. Note that the power spectrum depends only on the magnitude of the 
wavenumber, k = |k|, due to the assumptions of homogeneity and isotropy. 
Because 6, (a) x D(a) by definition, the dependence of the linear matter power 
spectrum on the linear growth is simple enough, 


P(k, a) x D(a)? x [ag(a)]?, (13) 


or, entirely equivalently, P(k,z) x D(z)*. Therefore, a simple observation 
of how galaxy clustering scales with redshift is sensitive to the growth of 
structure. This is, however, only one way to probe the growth, as cosmologi- 
cal observations are sensitive to its different aspects, only some of which are 
captured in Eq. (13). We now cover this in more detail. 


3 Connection to observations 


There are several ways in which one can measure the growth of struc- 
ture through cosmic time. Measuring galaxy clustering in redshift effectively 
measures P(k,z) in (typically) several redshift bins, thus probing growth. 
Additionally, one can map out the cosmic shear pattern, measuring correla- 
tions of shapes of distant galaxies which is sensitive to the distribution of 
mass along the line of sight and thus the growth of structure. Cross-correlating 
galaxy clustering and cosmic-shear signal is another method that can provide 
additional information on growth. A rather different approach entails mea- 
suring the space density of clusters of galaxies; their density as a function of 
mass and redshift — the mass function — depends on the linear growth func- 
tion D(z). Finally, one can measure the correlation of galaxy velocities as it 
too is sensitive to growth. We now discuss these cosmological probes, and how 
growth enters them. 


3.1 Galaxy clustering 


Because galaxy clustering is one of the oldest and best established cosmological 
probes, its measurements provide the principal way to isolate and constrain 
the growth of cosmic structure. Measurements of galaxy clustering constrain 
the power spectrum P(k, z) or, equivalently, the two-point correlation function 
&(r, z), over a range of scales and in several redshift bins. It is the dependence 
on redshift z that informs us about the temporal growth of cosmic structure. 

We now quantify the dependence of clustering on growth. The linear (|6| < 
1) dimensionless matter power spectrum A?(k,z) can be expressed in terms 
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ko P(k, z) 
Qn? 


So a EE) ep 
295 Q2; \kpiv Ho 1+z ? 


A?’ (k, z) = 


9 


where we expressed the temporal dependence in terms of redshift z. Here As 
is the normalization of the power spectrum (for the fiducial cosmology, As ~ 
2.1 x 107°), ns is the spectral index, kpiv is the “pivot” around? which A? (k) 
varies as a power law in wavenumber, the combination ag(a) = g(z)/(1+z) « 
D(z) determines the linear growth of perturbations, and T(k) is the linear 
transfer function which mainly encodes the change in the shape of the power 
spectrum around the scale corresponding to horizon size at matter-radiation 


equality. 


Dependence of the linear matter power spectrum on the (linear) growth 
in Eq. (14) is simple enough, P(k,z) « g(z)? x D(z)”. However, in practice 


there are two complications: 


e First, we can measure the power spectrum of galaxies and not dark-matter 


particles. Galaxy power spectrum is given in terms of the matter power 
spectrum P(k, z) as 

Pyalks2) = b?(k, z)P(k, z) (15) 
where b(k,z) is galaxy bias which depends on the wavenumber in a way 
that is difficult to predict theoretically. Further, galaxy bias depends on 
the galaxy type and, worse, on the galaxy formation history (often termed 
“assembly bias”), and thus typically needs to be measured directly from 
the data. Even on linear scales, where galaxy bias is expected to be scale- 
independent (i.e. constant in wavenumber k), its time dependence is a priori 
unknown. This time dependence of the bias is unfortunately degenerate with 
that of the growth of structure, and there is no way to break this degen- 
eracy with galaxy clustering alone. However, as we discuss in Sec. 3.2, one 
can use weak gravitational lensing (which altogether avoids galaxy bias), 
as well as galaxy-galaxy lensing (proportional to bias, rather than bias 
squared), to break this degeneracy and isolate the growth of structure. 
Alternatively, simultaneous measurements of the power spectrum and other 
statistics of galaxy clustering (e.g. the three-point correlation function, or 
the bispectrum) can help separately constrain galaxy bias and the growth 
of structure. 
The second complication is that typical clustering measurements are made 
partly in the quasi-linear and non-linear regime, corresponding roughly to 
scales r < 10h~'Mpe (or k = 0.1 h Mpc™t) at z = 0. Nonlinear corrections 


J: 


? Typically tajen to be kpiy = 0.05 Mpc ~, which is close to the wavenumber at which the 


primordial power is best constrained. 


Springer Nature 2021 ETRX template 


10 Growth of Structure 


to the power spectrum also depend on scale, but in a way that is not theo- 
retically tractable and must be calibrated with N-body simulations. Modern 
suites of simulations which predict the matter power spectrum at low red- 
shift do enable such a calibration, but the reliability of these predictions 
rapidly worsens as one goes to smaller scales. 


Traditionally, cosmological constraints from galaxy clustering have been 
quoted in terms of constraints on the amplitude of mass fluctuations o(z, R). 
This quantity effectively averages the matter power spectrum evaluated at 
some redshift z over a spherical region of comoving radius R. Mathematically, 


o?(z, R) = 7 A2(k, z) (EED) ang (16) 


where A?(z, k) is the dimensionless power spectrum from Eq. (14). The galaxy- 
clustering constraints are typically converted to those on the amplitude of mass 
fluctuations evaluated at z = 0 and R = 8h~!Mpc — hence we often make 
use of the quantity og = o(z = 0, R =8h~!Mpc). 

Galaxy clustering is a mature cosmological probe, as its first measurements 
date to 1960s and 70s. It is also reasonably easily accessible — one only needs 
to measure the galaxy positions in order to calculate the correlation function. 
Ideally galaxy redshifts are available, in which case full 3D power spectrum can 
be computed, otherwise one needs to determine the redshifts approximately 
using photometric information from galaxies (for a review of photometric red- 
shifts, see [53]). However, the comparison of galaxy clustering with theoretical 
predictions is seriously challenged by the presence of galaxy bias, as we dis- 
cussed above. This is where weak gravitational lensing comes in as a very 
powerful complement; we discuss this next. 


3.2 Cosmic shear 


Weak gravitational lensing — subtle distortions in galaxy shapes due to the 
intervening large-scale structure — is a very powerful probe of the growth of 
cosmic structure. Weak lensing is sensitive to the presence and distribution of 
mass along the line of sight, between the observers and source galaxies. This 
method, proposed in the 1960s and first detected in the year 2000, is now 
a standard-bearer for the probes of large-scale structure. Weak lensing also 
goes under the name cosmic shear as one statistically measures the amount 
of shearing of the observed shape of galaxies caused by photon deflections, 
which in turn informs about the projected mass along the line of sight. See 
Refs. [55, 56] for reviews of weak lensing / cosmic shear. 

The principal feature of cosmic shear is the absence of galaxy bias. Even 
though measuring galaxy shapes (for cosmic shear) is more challenging than 
measuring galaxy positions (for galaxy clustering), this absence of possible 
degeneracies between the bias and the growth of structure make cosmic shear 
a premier cosmological probe. 
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Fig. 2 Tomographic shear measurements of the 2-point correlation function of shear in real 
space, £7 (0) from the Year-3 data analysis of the Dark Energy Survey. The pair of numbers 
in € (written in the top left corner of each panel) refers to the pair of redshift bins from 
which the measurement is made (e.g. “2,4” means spatial correlation of galaxies in bin 2 to 


those in bin 4). The plots above the diagonal show oI (8) as a function of 8, while those 


below show ed (6). The curvy solid line shows the best-fit ACDM model. The grey regions 
show scales not used in the final parameter constraints due to concerns about modeling of 
the nonlinear clustering at these smaller scales. Adapted from DES Y3 shear analysis [54]. 


A key quantity in cosmic shear is the convergence «K, which is defined at 
every point on the sky and is proportional to the projected matter density 
between the observer and the source galaxy. Specifically, convergence is given 
by 


p= es f V?ð ds, (17) 
0 


where dz, dg, and dzs are respectively the distance from the observer to the 
lens and to the source, and the distance between the lens and the source. 
Further, ® is the three-dimensional gravitational potential that gets integrated 
along the line of sight in the above equation. A closely related quantity is shear 
y which, along with the convergence, makes up elements of a 2x2 matrix whose 
components define the distortion of an image at any point on the sky. 
Transforming the convergence to harmonic space and assuming statistical 
isotropy, one obtains the convergence power spectrum P*“", defined as (e.g. 
[57]) 
(kemkem) = Seer dmm! P (2), (18) 
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where the multipole £ correspond to an angle 0 ~ 180°/¢. The convergence 
power spectrum is identical to the shear power spectrum in the limit of weak 
distortions (i.e. in the weak-lensing limit); P77(@) ~ P**(@). For a Gaus- 
sian field, either PY or P** would contain all information. Beacuse the weak 
lensing field is nongaussian on small scales, higher-order correlations contain 
additional information and may be useful to exploit. 

The convergence (or shear) power spectrum is related to theory in a 
straightforward way: it is a projection along the line of sight of the three- 
dimensional matter power spectrum P(k). Here we also assume weak lensing 
tomography [58] — slicing of the shear signal in redshift bins — which enables 
extraction of additional information from the weak lensing shear, as it makes 
use of the radial information which is crucial for growth. Consider correlating 
shears in some redshift bin į to those in redshift bin j. Tomographic cross- 
power spectrum for these two redshift bins, at a given multipole £, can be 
related to theory by starting with Eqs. (17) and (18); the result is 


rosryo- [o ae (6). 


Here the weights W; are given by 


Wilx) = 5 MH aio (1 +2) (20) 
where a ( ) 
i = mi = 

ax) =r) | dmi S, (21) 


and n; is the normalized (f n(z)dz = 1) comoving density of galaxies if the 
coordinate distance x, falls in the distance range bounded by the ith redshift 
bin, and zero otherwise. Further, r(z) is the comoving distance and H(z) is the 
Hubble parameter. Note that the weights W; are purely geometrical, and do 
not contain any growth information. The key takeaway from Eq. (19) is that 
weak-lensing convergence power spectrum is proportional to the matter power 
spectrum P(k), and hence (on linear scales) to the growth squared, D(z)?, but 
without the complicating presence of galaxy bias. 

Weak lensing therefore offers particularly good prospects for constraining 
the growth of structure. With N tomographic bins, one has a total of N(N + 
1)/2 tomographic bin pairs, each of which contains cosmological information. 
For typical surveys N ~ 5-10, enabling in principle the constraints on the 
growth of structure D(z) with relatively fine temporal resolution. 

There are also important systematic errors that enter weak-lensing mea- 
surements. Notably, shear measurements are in general challenging, and are 
impacted by smearing of galaxy shapes by atmospheric distortions or other 
instrumental artifacts. Intrinsic alignments between source galaxies complicate 
theoretical predictions and need to be explicitly modeled. Finally, photomet- 
ric redshifts of source (and lens) galaxies have complex statistical properties 
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that need to be corrected for where possible, then explicitly modeled in the 
analysis. For a review of these and other issues, see [59]. 


3.3 Galaxy-shear cross-correlations and 3 x 2 analysis 


In addition to measuring shapes of distant galaxies across the sky and cor- 
relating them, there are other ways to leverage weak lensing observations to 
learn about the growth of structure in the universe. One option is to mea- 
sure the correlation of the shears of background galaxies with the positions 
of the foreground galaxies. This correlation is known under the name galasy- 
galaxy lensing, and effectively measures the lensing efficiency of foreground 
galaxies. The method should perhaps be called galaxy-galaxies lensing, as it 
correlates the position of one foreground galaxy with the shapes of a number 
of background galaxies that are near it on the sky. Another better name is 
galaxy-shear cross-correlation, as galaxy positions are correlated with (other 
galaxies’) shear. 

Modern analyses combine the galaxy clustering (galaxy-galaxy correla- 
tions), cosmic shear (shear-shear correlations), and galaxy-galaxy lensing 
(galaxy-shear correlations) into one unified analysis. This combination goes 
under the informal name of “3 x 2” analysis, named for three two-point cor- 
relation functions. Schematically, 3 x 2 observations can be organized into a 
matrix 


3 x 2 analysis ~ oe , (22) 
IY YY 


where g stands for galaxy position, and y for cosmic shear. If the actual obser- 
vations were made in a single redshift bin and on a single spatial scale, then 
they would indeed correspond to the data vector of three elements, [gg, gy, y7]. 
Because the observations span a number of redshifts and scales, the resulting 
data vector can be much longer. For example, in the DES Y3 analysis [60], 
the full data vector contains about 1000 elements, corresponding to measure- 
ments of gg, gy, or yy, evaluated at a number of spatial scales k and in several 
redshift bins. About half of these measurements were made on linear or quasi- 
linear scales where theoretical modeling is accurate and were used to constrain 
cosmological parameters. 

Because a 3 x 2 analysis efficiently combines observations of galaxy clus- 
tering and gravitational lensing, it also serves as a major probe of the growth 
of structure. The constraints on growth from a 3 x 2-type analysis are not 
trivial to isolate; we will discuss in Sec. 4.2 how this can be done. Instead, we 
typically quote constraints on the combination of the amplitude of mass fluc- 
tuations and matter density that is well constrained by galaxy clustering and 


weak lensing 
Qu 
Sg = —. 23 
8 = e853 (23) 
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Constraints on Sg are then roughly interpreted as those on the growth of cosmic 
structure. In Sec. 3.7 below, we show a compilation of current constraints on 
Sg from current data. 


3.4 Counts of galaxy clusters 


The abundance of clusters of galaxies as a function of their mass and redshift 
is very sensitive to cosmic growth. Galaxy clusters are the observed baryonic 
components of dark-matter halos, so the abundance of clusters can often be 
quantified via our theoretical understanding of the abundance of halos. The 
halo abundance formalism can be studied from first principles using the excur- 
sion set formalism pioneered by Bond et al [61] and reviewed by Zentner [62]. 
In this formalism, formation of a halo occurs when random walk of overden- 
sity 6, as a function of radius of a sphere over which it is evaluated, crosses 
some critical threshold. This threshold is quoted in terms of the peak height 
v = e/o, where ĝe ~ 1.686 is the critical overdensity for collapse and ø is the 
amplitude of mass fluctuations defined in Eq. (16). Because halo formation 
is exponentially sensitive to the peak height (in the simple spherical-top-hat- 
overdensity model), and v x 1/a(z) x 1/D(z), the sensitivity of the halo 
abundance to the growth of structure is very strong. 

Because galaxy clusters are the most massive collapsed objects in the uni- 
verse, they are described by the density field in the not-overly nonlinear regime. 
As a consequence, their abundance can be modeled analytically, using argu- 
ments first laid out by Press and Schechter [63] that preceded those of the 
more general (and aforementioned) excursion-set formalism, and with a rather 
transparent dependency on the cosmological parameters [64, 65]. The quantity 
that describes the abundance of galaxy clusters is the mass function, which 
is the space density at arbitrary redshift and mass, n(z, M), and is typically 
reported as the number per unit log interval in mass, dn/dln M. The number 
of clusters with mass above Mmin and below redshift Zmax in some volume of 
solid angle Qsky is then 


Zmax co 2 
N(z < zmax, M > Mmin) = Qsky f dz / a?) na, G0 
0 Mmi 


an dn M H(z) 


where r?(z)/H(z) = dV/(dQdz) is the volume element, with r the comoving 
distance and H the Hubble parameter. In practice, the number of clusters 
is evaluated by connecting it explicitly to directly observable quantities (as 
discussed near the end of this subsection), but the mass function retains its 
central role in connecting observations with theory, and its strong dependence 
on the growth of structure remains useful. 

While the cutting-edge approach for modeling the mass function is to 
calibrate it directly on simulations and interpolate in space of cosmological 
parameters using a so-called emulator [67-69], here we stick with the semi- 
analytical mass function for illustrative purposes. We adopt the Tinker [66] 
mass function, which has been shown to be accurate to ~ 5% over a respectably 


Springer Nature 2021 PTẸX template 


Growth of Structure 15 

= solid: Qy = 0.30 5 f 
> wl ee dashed: Q y =0.25 ae 
a, ~~ ~ —— QNyu=0.25 
a Sam pagan ge 
z So === 28" = 0.30, 04" =0.25 
~ 
2 © E 
= z3 
T = 
x S Vary growth 
= A 2T and volume 
3 S — 
= ZIE O aie 
> Vary growth alone 
3 N, 

1010 1 1 14 poiitiiit Lios P O O S ee Se 

0.0 . ` ; 4 $ 3.0 8o 0.2 0.4 0.6 0.8 1.0 
z z 
Fig. 3 Sensitivity of galaxy cluster counts on growth. The left panel shows the Tinker 


et al. [66] mass function, as a function of redshift, for M = 1018, 1014, and 1015 Mo. The 
right panel shows the number of clusters expected above M = 10!4Mo out to redshift z, 
assuming a survey covering 5000 sq. deg (note that the y-scale shows counts in units of 104). 
The right panel shows that the sensitivity of number counts on cosmology comes largely 
through growth, since changing the growth alone (and keeping geometrical terms unchanged) 
accounts for most of the difference when moving to a new cosmological model — here, a 
different value of Qm. 


wide range of mass and redshift. The Tinker mass function is defined as dM 


and dz as 
dn 


dinM 
aaler fem 


where pmo is the matter density today, ø is the amplitude of mass fluctuations, 
and a,b, and c are redshift-dependent coefficients calibrated from simula- 
tions and reported in Ref. [66]. Moreover, M refers to halo mass specifically 
defined as the mass in a spherical region with 200 times the mean mass den- 
sity in the universe. Here, the growth function dependence enters through the 
amplitude of mass fluctuations o(M, z) œ D(z) which is found in the exponen- 
tial, as predicted by the Press-Schechter (and, more generally, excursion-set) 
arguments. 

In Figure 3 we show the sensitivity of the cluster counts on the growth 
of structure. The left panel shows the Tinker mass function vs. redshift for 
three values of mass (M = 10!8,10!*, and 10! Mọ), and for two values of 
matter density Qw. It shows strong dependence of the expected number den- 
sity of clusters on the matter density, which in turn determines the growth 
of structure. The right panel shows the number of clusters expected above 
M =10'*Mo out to redshift z, assuming a survey covering 5000 sq. deg (note 
that the y-scale shows counts in units of 104). We observe that the sensitiv- 
ity of number counts on cosmology comes largely through growth, rather than 
the geometrical volume factor. This explicitly illustrates the fact that cluster 
counts are very sensitive probes of growth [65, 70-76] 

Because cluster masses are not directly measured, modern measurements 
of the abundance of clusters are typically compared to theory not in terms of 


pmo dino! 
ea ane 


(25) 
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their masses, but rather other intermediate, more readily observable quanti- 
ties called mass “proxies”. One such mass proxy is “richness”, defined as the 
number of galaxies that reside in a cluster. Other mass proxies include X-ray 
flux or weak-lensing signal measured toward clusters; all of these proxies cor- 
relate with cluster mass. Conversion from noisy measurements of mass proxies 
to actual cluster masses introduces both statistical and systematic errors, and 
controlling and quantifying these errors — especially the systematics — is the 
principal challenge for cluster cosmology. If the systematics can be controlled 
and understood, then the prospects for constraining the growth of structure 
via cluster abundance with ongoing or upcoming wide-field cosmological sur- 
veys such as eROSITA, LSST, and Euclid and Nancy Roman Space Telescopes 
are very good. 


3.5 Cosmic velocities 


Cosmic velocities are also sensitive to the growth of structure. In linear theory, 
the velocity v is directly related to overdensity 6 via the continuity equation 
which reads (see e.g. Ref. [77]) 
oy, Say i 

v= aoa af He, (26) 
where the prime denotes a derivative with respect to conformal time 7, and 
the second equality follows because ’ = d/dn = a(d/dt) = aH(d/dlna). Here 
we have defined the growth rate f as 


_ dinD 


f= dlna` 


(27) 


Then the velocity power spectrum can be related to the matter power spectrum 
in linear theory via 


Pos(k,a) = [00] P(k, a). (28) 


Note a key feature: the velocity power spectrum not only scales as the matter 
power spectrum (and hence the usual growth term D(a) squared), but is fur- 
ther proportional to the square of the growth rate f(a). Because the growth 
rate is sensitive to both the standard cosmological parameters (Qm and w, for 
example) and modified gravity, this latter dependence makes the cosmic veloc- 
ities particularly well suited to cosmological tests that rely on the growth of 
structure [78]. 

Probably the most prominent method for utilizing the velocities are the 
redshift-space distortions (RSD). The RSD take place because redshift of 
galaxies is affected by their gravitational infall into nearby large-scale struc- 
tures, as well as by the galaxies’ own peculiar velocities. Because we typically 
measure the galaxies’ redshifts in order to get their radial location, clustering 
measurements in redshift space are subject to RSD. There are two RSD effects 
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Fig. 4 Illustration of the redshift-space distortions. There are two principal effects: 1) On 
large scales, velocity flows into large overdensities, “squishing” the appearance of the object 
along the line of sight; this is the Kaiser effect; 2) On smaller scales, random motions 
introduce apparent elongation along the line of sight; this is the somewhat hyperbolically 
called “fingers of god” effect. Adopted from [57]. 


— the so-called Kaiser effect, and fingers of god; see Fig. 4. To lowest order, 
the redshift-space matter power spectrum P(k)*) is related to the isotropic 
matter power spectrum P(k) via (e.g. [78])® 


P(k, 2)® = [b+ fu)” F(k?o2p?)P(k, 2), (29) 


where the superscript ‘*) indicates that we are referring to clustering in redshift 
space. Here y is the cosine of the angle between the line-of-sight direction and 
wavenumber k, b = b(k,z) is the galaxy bias, f = f(z) is the growth rate, 
and P(k,z) is the isotropic part of the matter power spectrum. Finally, the 
function F(k?02?) ~ 1/(1 + k?o2u?) models the suppression of the redshift- 
space power spectrum at high k in terms of the velocity dispersion o», which 
is a free parameter that generally depends on halo mass and redshift. 

The RSD effects enter on mildly nonlinear scales (k > 0.1 h Mpc~!), where 
they can be readily measured in spectroscopic surveys where accurate galaxy 
redshifts are available. Because P(k)‘*) has terms that go as y°, i, and u4 


3A slightly more accurate formulation replaces the density power spectrum P(k,z) = Pss on 
the right-hand side of Eq. (29) with a combination of Ps5, the velocity power spectrum Psv, and 
the cross-power P5», as [79] 


P(k, 2) = [b? Pss Jofa Psy + Pur Pao] F(k?02 7). 


Springer Nature 2021 TeX template 


18 Growth of Structure 


(in linear theory), we can respectively measure the monopole, quadrupole, and 
hexadecapole of redshift-space galaxy clustering. 

An alternative to utilizing the RSD to probe growth is to measure the 
velocities directly. Redshift to a galaxy is affected by its velocity, so accurate 
measurements of the former can constrain the latter, and hence cosmic growth. 
Velocity of a galaxy at a physical distance x is 


= L (ar) = år + at = HX + Vpec; (30) 


x 
where r is the comoving distance, H is the Hubble parameter, and Vpec is the 
component of the peculiar velocity of the galaxy parallel to the line of sight. 
Specializing in z < 1 (and ignoring 2nd-order terms in redshift), we get the 
relation between the observed redshift zp, and true redshift z as 


CZobs = CZ + Upec: (31) 


The effect of peculiar velocities is non-negligible only for very nearby galaxies; 
consider that for a galaxy at z ~ 0.03, cz ~ 10,000kms~', while the typical 
peculiar velocity is Upec ~ 200kms™!. 

To determine the peculiar velocity, we need the observed redshift and 
distance. Galaxy redshift can be determined by taking a sufficiently high- 
resolution spectrum of the galaxy, and is in principle straightforward. Measur- 
ing the distance to a galaxy is far more challenging. To get the distance, one 
may use an empirical fundamental-plane relation that relates an object’s dis- 
tance to its surface brightness. Alternatively, one can use type Ia supernovae, 
standard candles whose distances can be determined by measuring their fluxes 
and adopting the fact that these objects have nearly equal luminosities. In sur- 
veys thus far, the fundamental-plane measurements are available for ~10,000 
galaxies and give distances accurate to 20 — 30% per object. Type Ia super- 
novae are more accurate (~10% per object), but available for fewer galaxies 
(typically < 500) in extant peculiar-velocity surveys. Either method allows 
a noisy measurement of individual galaxies’ peculiar velocities; their correla- 
tions are sensitive to growth as indicated in Eq. (28). Recent observational 
efforts, combined with modern distance calibrations, have been greatly increas- 
ing the number of galaxies and supernovae available for peculiar-velocity tests 
of cosmology [80, 81]. 

For practical reasons having to do with patchy sky coverage and sparsity 
of source galaxies, one often measures the (configuration-space) correlation 
function of galaxy velocities. The correlation function Cy, is proportional to 
the velocity power spectrum. Since the velocity power spectrum is proportional 
to the quantity fos = f(z)o8(z) 


Polk, z) x f? (2) P(k, 2) = P (z)o8(z)P(k, z = 0), (32) 
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so is the velocity correlation function 
Cov x Puy x (fos). (33) 


This relation explains why the results from the analyses of cosmic velocities 
are often framed as constraints on the quantity fog, evaluated in redshift bins. 
Note however that, while the strongest dependence of velocity correlations is 
on the quantity fog, most of the standard cosmological parameters also enter, 
mainly via their impact on P(k) that is not captured by og. This complicates 
attempts to isolate the measurements on fog, as these other parameters (e.g. 
the physical baryon density Nph? or the spectral index n,) need to be given 
priors from other data, or else be simultaneously constrained by the velocities. 

Observing the peculiar velocities of nearby galaxies or type Ia supernovae 
that they host is not the only way to probe the galaxy velocity field. One 
useful probe in this regard is the kinetic SZ effect (kSZ). The kSZ is caused 
by galaxy clusters’ peculiar velocities in the CMB rest frame. A cluster with 
peculiar velocity Upec, along the line of sight will incur a density shift of the 
CMB temperature in the direction of the cluster of 


(=) ~ —Te pec, | ; (34) 
T J ksz ç 


where Te is the optical depth for electron scattering in the cluster. The kSZ 
modifies the blackbody spectrum of an object by the temperature shift given 
in Eq. (34). For typical values Te ~ 0.01 and Upec ~ 500kms~', the kSZ 
fractional temperature shift is on the order of the primordial temperature 
anisotropy (~ 1075) and thus quite small. The kSZ effect was first detected 
about a decade ago [82], and its better mapping will allow us to probe the 
velocity field of galaxy clusters, and hence the growth of cosmic structure. 


3.6 Role of the CMB 


Being a snapshot of the universe at z ~ 1000, the cosmic microwave back- 
ground (CMB) naively does not have much information about the growth of 
structure. Specifically, we are most interested in the growth at relatively low 
redshift, z < 2, where growth slowly transitions from being robust in the dark- 
matter dominated regime, to being strongly suppressed as dark energy starts 
to dominate*, and we certainly do not expect that the CMB be informative 
at such a low redshift. While technically correct, these expectations need to 
be modulated with the fact that the CMB is nevertheless extremely impor- 
tant in pinning down the initial amplitude of the power spectrum (that is, the 
parameter A, in Eq. (14)). This constraint, along with the present-day clus- 
tering amplitude (measured as og or Sg), allows a significant improvement in 
the precision of constraints on the growth of structure. 


4In the standard ACDM model with 30% matter and 70% dark energy, the energy densities in 
these two components are equal at z ~ 0.33. 
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Fig. 5 Constraints on fog = f(z)og(z). The curves show predictions from ACDM with two 
alternate values of the matter density — the fiducial Qm = 0.3, and also Qm = 0.25. The 
errors show the various extant measurements, as follows. The black error bars denote con- 
straints from peculiar velocities (or their combination with clustering) at z œ 0; these black 
points are mutually offset horizontally for clarity. From left to right, these peculiar-velocity 
constraints come from: type Ia supernovae from the 2M++ compilation [87], supernova dis- 
tances from Supercal compilation and Tully-Fisher distances from 6dFGRS [88], combination 
of galaxy velocities and galaxy clustering applied to 6dFGRS [89], and combination of galaxy 
velocities from 6dFGRS and SDSS with velocity and density predictions from 2M+-+ galaxy 
survey [90]. The other constraints, at higher redshift, come from redshift-space distortion 
modeling of galaxy samples; see text for details. 


CMB contributes in other ways. Notably, CMB lensing — measured as 
subtle displacements in the distribution of hot and cold spots on arcminute 
scales — is sensitive to the growth of structure. CMB lensing is described by 
the power spectrum of the deflection signal of photons as they cross large-scale 
structure traveling toward our detectors. Mathematically, the CMB lensing 
power spectrum looks similar to the weak-lensing shear power spectrum in 
Eq. (19); the main difference is that kernel of the shear power spectrum peaks 
at z ~ 0.5, while that of CMB lensing peaks at z ~ 3. Hence, CMB lensing in 
principle helps probe the growth of structure at a redshift higher than galaxy or 
weak-lensing surveys. The signal-to-noise of CMB lensing is currently limited, 
but future constraints from surveys like Simons Observatory and CMB-S4 
may provide interesting constraints on the growth of structure in their own 
right. Another promising direction is to cross-correlate the CMB map with a 
galaxy map; the contribution of CMB lensing should lead to a sufficiently high 
signal-to-noise to probe the growth of structure [58, 83-86]. 
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3.7 Cosmological constraints on fog and Ss 


Having reviewed the principal cosmological probes of growth, we now summa- 
rize current constraints. While the temporal evolution growth enters many of 
the observable quantities in cosmology, isolating it is in general not straight- 
forward; we discuss such strategies in Section 4 below. Instead, the majority 
of the cosmological analyses assume that growth is modeled according to the 
assumed cosmological model (e.g. ACDM and its variants, or else modified 
gravity), and report the constraints on the parameters of that model, say Qm 
and og. In order to present results in a more “model-independent” way, these 
surveys report constraints on derived cosmological parameters that are fairly 
directly related to what is being measured, yet contain information that in 
large part comes from growth. We now discuss these constraints. 

Perhaps the most direct constraint comes from peculiar-velocity and 
redshift-space-distortion measurements, both of which are sensitive to the 
quantity fog at some effective redshift that depends on the distribution of 
source galaxies (see Eqs. (28) and (33)). The quantity fog depends on both 
growth and other cosmological parameters, principally via how they enter the 
matter power spectrum. Despite some degeneracy with these other parame- 
ters, the dependence on growth is very strong because both f(z) and og(z) are 
directly related to the linear growth D(z); fog x (dln D/dlna)D. 

Constraints on fog from current measurements are shown in Fig. 5. The 
black error bars denote constraints from peculiar velocities (or their combina- 
tion with clustering) at z ~ 0; the points are mutually offset horizontally for 
clarity. From left to right, these peculiar-velocity constraints come from: type Ia 
supernovae from the 2M+-+ compilation [87], supernova distances from Super- 
cal compilation and Tully-Fisher distances from 6dFGRS [88], combination of 
galaxy velocities and galaxy clustering applied to 6dFGRS [89], and combi- 
nation of galaxy velocities from 6dFGRS and SDSS with velocity and density 
predictions from 2M++ galaxy survey [90]. The other constraints, at higher 
redshift, are mainly from redshift-space distortion modeling of galaxy samples. 
These RSD constraints come from: 6dFGRS [91], GAMA [92], WiggleZ [93], 
VIPERS [94], SDSS main galaxy sample (MGS) [95], FastSound on Subaru 
[96], and BOSS+eBOSS [97]. We see that the fog constraints are in a good 
agreement with the predictions of the currently favored ACDM model. Partic- 
ularly notable is the complementarity between the different probes, as peculiar 
velocity surveys measure motions of galaxies at distances r < 100 h~'Mpc, and 
thus very low redshift (z ~ 0.02), while the RSD measurements probe growth 
to galaxies and quasars at a much higher redshift, all the way up to z ~ 1.5 
for current surveys. 

In contrast to peculiar-velocity and RSD measurements, constraints from 
weak-lensing surveys and measurements of the broadband galaxy-clustering 
power spectrum are not very sensitive to the growth rate f or the combination 
fog, but rather on the overall amplitude of matter fluctuations. This amplitude 
is larger if either matter density Qa, or the amplitude of mass fluctuations og 
increases. Therefore, the constraints typically look like a banana-shaped region 
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Fig. 6 Constraints on Sg = og,/Qy,/0.3 from various analyses. Going top to bottom, 
the first three measurements come from the CMB analysis: temperature-only analysis from 
Planck ([98], also shown as a vertical shaded region), temperature + polarization + lensing 
analyses from Planck [98], and temperature + polarization information from ACT combined 
with large-scale temperature from WMAP [99]. All of the other constraints involve some 
combination of large-scale structure probes, see text for a detailed description. All constraints 
assume the fiducial ACDM model. The error bars show 68% credible intervals, which are 
asymmetric around the mean in some cases. 


in the Qm — og plane, indicating their mutual anti-correlation. To decouple 
these two parameters, the constraints are often reported on their combina- 
tion Sg = og(Qu/0.3)°° that is very well constrained. Note that og and Sg 
explicitly refer to the value of these quantities at redshift zero, to which the 
weak-lensing and galaxy-clustering are extrapolated by convention, in contrast 
to peculiar-velocity and RSD constraints which are traditionally quoted as 
fos = f(z)os(z). 

Figure 6 shows the constraints on Sg. Going top to bottom, the first 
three measurements come from CMB analyses: temperature-only analy- 
sis from Planck ([98]; also shown as a vertical shaded region), tem- 
perature+polarization+lensing analysis from Planck [98], and tempera- 
ture+polarization information from ACT combined with large-scale tempera- 
ture from WMAP [99]. Going further down, we show the constraints of select 
large-scale structure combinations as follows. The black error bars show con- 
straints from the fiducial 3x2 analysis from DES Y3 [22], DES Y3 with cosmic 
shear alone [54, 100], and DES 3x2 with the addition of Planck lensing infor- 
mation [101]. The purple error bars show constraints from KiDS-1000 [21] and 
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Fig. 7 Constraints on the time-dependent amplitude of mass fluctuation og(z) and com- 
parison with theory, adopted from DES Y3 extensions paper [22]. The blue points show 
the constraints from DES 3 x 3 analysis alone in each of the four lens redshift bins. The 
black points show constraints from the combination of DES 3 x 2 and external data (Planck, 
BAO, RSD, SN), along with the shear ratio statistic from DES (labeled SR). Lines and 
shaded bands show the means and 68% credible intervals inferred from ACDM posteriors 
corresponding respectively to DES alone (blue), and DES plus external data (black). 


KiDS-1000 combined with clustering from BOSS and galaxy-galaxy lensing in 
the overlap region of KiDS, BOSS and 2dFLens [102]. The green error bars 
show the constraints from harmonic-space [23] and real-space [103] analyses of 
Hyper Suprime-Cam (HSC) year-1 data. The red error bars show constraints 
from the combined BAO and RSD information from BOSS+eBOSS [97], from 
the emulator-based full-shape analysis from BOSS DR 12 [104], full-shape 
BOSS DR12 analysis combined with the BOSS DR12 bispectrum monopole 
[105], and full-shape BOSS DR12 along with its cross-correlations with Planck 
lensing [106]. The brown error bars show combination of unWISE galaxy clus- 
tering and Planck lensing [107], and clustering of luminous red galaxies from 
DESI imaging survey combined with Planck lensing [108]. The turquoise error 
bars show constraints from the abundance of galaxy clusters as measured by 
South Pole Telescope [109] and DES Year-1 data [110]. All constraints assume 
the fiducial ACDM model. 

The most apparent trend in Fig. 6 is the so-called? “S8 tension”, which 
reflects the fact that CMB measurements show a higher amplitude of mass fluc- 
tuations than lensing surveys. Planck’s temperature, polarization, and lensing 
information combined indicate Sg = 0.832+0.013, and the combination of ACT 
and WMAP give a consistent result (Sg = 0.840 + 0.040). Lensing surveys, 
on the other hand, typically show a lower value of this parameter, exempli- 
fied by DES Y3 constraint Sg = 0.775 + 0.017 and a very similar constraint 
from KiDS (Sg = 0.75970')3+). Similar trends are seen with cluster abundance 
constraints. Future data will sharply improve the constraints on the growth. 

Constraints on og or Sg are typically projected to, and reported at, red- 
shift zero. One can do better however, and constrain the temporal evolution of 
the (linear) amplitude of mass fluctuations, with the understanding that these 


5 . > . . . . . 
? Also referred to as the “og tension”, as a similar trend is seen in constraints on og. 
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amplitude parameters are directly proportional to the linear growth function 
D(z). Such constraints are shown in Figure 7, adopted from DES Y3 extended- 
model paper [22]. The blue points show the constraints from DES 3 x 3 analysis 
alone in each of the four lens redshift bins. The black points show constraints 
from the combination of DES and external data: CMB, baryon acoustic oscil- 
lations (BAO), redshift-space distortions (RSD), and type Ia supernovae (SN). 
Lines and shaded bands show the mean and 68% credible interval inferred 
from ACDM posteriors corresponding respectively to DES alone, and DES 
plus external data. Clearly, og(z) measurements are already accurate, espe- 
cially when various datasets are combined, and sharply test (and, thus far, are 
in agreement with) the standard ACDM cosmological model. 


4 Consistency tests with growth 


As mentioned in the Introduction, growth of structure owes its outsize impor- 
tance chiefly to the fact that it is a powerful discriminator between models of 
dark energy, and that it is sensitive to the presence of modified gravity (e.g. 
[111, 112]). In Section 3 we discussed how individual probes are sensitive to 
growth, and in particular (in Sec. 3.7) how the contraints derived from growth 
are reported in the measurements of fog(z) and Sg. We now discuss more 
ambitious and direct uses of the information stored in the growth of structure, 
specifically how to utilize it to test classes of cosmological models. 


4.1 Constraints on parameterized growth 


While it is possible to isolate the constraints coming specifically from the 
growth of cosmic structure (as we discuss below in Sec. 4.2), this is typically not 
straightforward to implement. Moreover, such relatively general extractions of 
the growth information necessarily impose additional parameters that weaken 
the cosmological constraints. This is where a simple parameterization of growth 
becomes extremely useful. A successful such parameterization would have the 
ability to describe growth in ACDM and wCDM (where the equation of state 
of dark energy w is a free but constant parameter), as well as model departures 
expected in modified-gravity scenarios. 

By far the most impactful parameterization of this kind is provided by 
“growth index” parameter? y. This description of growth introduces a single, 
constant parameter defined in a phenomenological fit to the growth rate [113] 


f(a) = Nmulay. (35) 


In other words, the linear growth factor is approximated by D(a) = 
efo dina’ [Qu (a)"], 


It has been known for a long time that the formula in Eq. (35) fits the 
linear growth very well for y ~ 0.55 (e.g. [114, 115]). What is new here is 
that promoting y to a free parameter enables describing growth in models far 


Not to be confused with shear which we discussed in Sec. 3.2. 
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beyond the standard ACDM. Specifically, it has been shown that dynamical 
dark energy models where the equation of state ratio is parametrized as w(a) = 
Wo +Wa(1—a) are well fit by the growth index as long as it takes the value [113] 


y = 0.55 + 0.02[1 + w(z = DJ, (36) 


where w(z = 1) is the dark-energy equation of state evaluated at redshift one 
in this class of models. The form in Eq. (36) fits the exact linear rate to better 
than 0.3% when wo + wa < —0.1 [116]. A broad range of modified-gravity 
models, including time-varying gravity, DGP braneworld gravity [117], and 
some scalar-tensor gravity, are fit accurately with the growth index [118]. 

There are other variants of parametrized growth. A more direct, and even 
less model-dependent, approach is to model D(z) as a free function in redshift 
and interpolate it with principal components [58], splines [86], or else with 
piecewise-constant values as done for og(z) in the DES Y3 extensions paper 
([22]; see our Fig. 7). As the measurements of the growth of structure become 
more precise with upcoming surveys, such ambitious approaches will begin to 
return very interesting constraints and consistency-test results. 


4.2 Comparing growth with geometry 


Comparing measurements of geometric quantities to those describing the 
growth of structure is a particularly promising stress-test of the cosmological 
model [116, 119-123]. For example, given constraints on the initial conditions 
(power spectrum shape and amplitude), very precise distance measurements 
from e.g. type Ia supernovae (SNe Ia) and baryon acoustic oscillations (BAO) 
predict the convergence power spectrum measured by weak lensing probes. 
Here, the weak lensing signal depends on the late-time growth of structure, 
which in turn is precisely determined by distance measurements. 

A straightforward way to compare geometry and growth, first proposed in 
the modern form and applied to early data by [124], and then further developed 
by [125-129], is to add additional cosmological parameters. Specifically, in the 
flat ACDM cosmological model one can take the matter density relative to the 
critical — which is normally encoded in one parameter, Qw — and duplicate 
it into two parameters: Qf" and 05°". One then modifies the theory code as 
follows: every geometrical term (for example, a distance formula) will depend 
on Qf; ", and every growth term (for example, in the growth-of-structure 
differential equation) will be fed 047°”. With such an implementation, the tight 
relations linking geometry to growth in the standard cosmological model will 
be explicitly decoupled. In a flat cosmological model, the energy density of dark 
energy is given by Qa = 1—Qy, so such a split of Qm is also automatically a 
split of Q4. One can extend such a geometry-growth split to models with more 
complicated dark-energy sectors. For example, in the analogous scenario of the 
flat wCDM model, the equation-of-state of dark energy, w, is also described 
by two parameters, w®°°™ and w®°’. Therefore, a cosmological analysis with 
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Fig. 8 Geometry-growth constraints from individual probes as well as the joint analysis. 


Here we extend the wCDM model by assuming two parameters for the matter density relative 


to critical, Q$% 7, which governs the gemetrical quantities, and QO”, which controls the 


growth of structure; we also split the equation of state into w8°°™ and w8'°”, and marginalize 
over both of these parameters (and all others, unsplit parameters) in this Figure. The dashed 


black line shows the consistency relation that is satisfied in the standard cosmological model, 


OF = Q7- Note the impressive degeneracy breaking when the cosmological probes are 


combined. Adopted from Ref. [129]. 


geometry-growth split is specified by 


e split ACDM: {080™, 08" {p,}} (37) 


e split wCDM : {osom genr. wee wW, {pi} (38) 


where {p;} are other, standard (and unsplit) cosmological parameters. 

Constraints on the geometry and growth parameters in the split wCDM 
model are shown in Fig. 8, adopted from [129]. Here we show the constraints 
on OF" and Qf, marginalized over all other parameters, including w®°°™ 
and wystow, Note that SN Ia are manifestly sensitive to geometry only, and 
not to growth. The constraints from individual probes are necessarily weak, 
as the split parameter space contains severe parameter degeneracies. When 
the probes are combined, however, parameter degeneracies are broken and the 
constraints are strong. 

Current constraints on geometry and growth indicate intriguing departures 
from the standard cosmological model (which predicts QGP" = OQ” and 
weeem = wW) with some combinations of datasets (125, 126, 129). These 
results will be well worth following up with forthcoming data. 


4.3 Falsifying classes of models 


Another very effective way to test the consistency of any given cosmological 
model is to: 
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e Start with constraints on the model that are provided by certain datasets 
or probes; 

e Compute what those parameter constraints imply for the larger space of 
observable quantities, for example cosmological functions D(z), or fog(z) 
at an arbitrary redshift, or else derived parameters such as Sg; and then 

e Test those predictions using new data or probes. 


Such an approach has implicitly been discussed for a long time (e.g. [116]), but 
was first clearly spelled out by Mortonson et al. [130] and applied to current 
data by Refs. [131-134]. 

Falsification procedure as described above has a few clear advantages. First, 
the simple structure of smooth dark-energy models implies that, given cur- 
rent constraints, these models lead to very precise predictions for the temporal 
evolution of distance and growth at all redshifts, even well beyond the range 
probed by data used to impose the predictions. Second, those predictions are 
quite accurate not only for ACDM and wCDM, but also for models with a 
much richer dark-energy sector [130]. This makes the predictions good tar- 
gets for falsifying whole classes of models, particularly with data that probe 
new ranges of redshift (or new spatial scales) than those used to impose the 
predictions. Third, the ability to separately test growth and geometry comes 
about naturally in these tests, as predictions for both quantities as a function 
of redshift can be made straightforwardly, essentially by running the standard 
theoretical calculation constrained to the range of cosmological parameters 
allowed by current data. 

We have already shown one example of predictions in this review; this is 
the Planck measurement of Sg, shown as the vertical band in Fig. 6. Note that 
the spirit of model predictions is at work here: even though Ss is not a core 
parameter that CMB constrains, its value can be easily predicted by the CMB 
anisotropy in the given cosmological model (here, ACDM). This measurement 
serves as a check against more direct Sg measurements from galaxy clustering 
and weak lensing. 

There are other ways to falsify models that are not quite as direct as the 
procedure described above. One such test is to compare the two gravitational 
potentials V and ®. Starting from the scalar-perturbed Friedmann-Robertson- 
Walker line element in the conformal Newtonian gauge 


ds? = a?(n) [(1 + 2W)dn? — (1 — 26)6,;daide,] , (39) 


the behavior of many cosmologically interesting modified-gravity theories can 
be fit with a free function of time and scale multiplying the Poisson equation, 
and another given by the ratio of ® and W. Specifically, one introduces two 
new functions of time u(a) and U(a), defined as 


k?’ V = —4nGa?(1 + p(a)) 6 
(40) 
k?(W + ©) = —87Ga?(1 + D(a))p6, 
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where ô is the comoving-gauge density perturbation. Here,  parametrizes the 
change in the lensing effect on massless particles for a given matter field, while 
u describes the change in the matter overdensity itself. Such a parameterization 
approximately describes a more complicated set of equations [135-144], but 
has been numerically verified on scales of cosmological interest [45, 47-49, 145]. 
Typically, one either assumes that u and ÈX are constant, or else that they 
scale with cosmic time under some parameterization [146, 147]. Finally, there 
exist other, closely related, parameterizations of the gravitational potentials, 
including the Eg parameterization [148] and parameterized post-Friedmannian 
framework [149]. 

Because X determines the predictions for lensing, weak lensing measure- 
ments are primarily sensitive to this parameter but also have some smaller 
degree of sensitivity to u via their tracing of the matter field. Conversely, galaxy 
clustering measurements depend only on yz and are insensitive to © [147]. Con- 
straints on (u, X), or similar quantities but with different names, are already 
very good [147, 150-152]. More ambitious analyses promote u and © to few- 
parameter functions of redshift [22, 98, 147, 153-155]. Comparison of (a) and 
(a) to their fiducial values of zero thus constitutes a test of modified gravity. 


5 Conclusions 


We have described how the growth of structure determines key aspects of 
observable quantities in cosmology, such as the spatial correlation of galaxies, 
the coherence of their peculiar velocities, the amplitude and spatial depen- 
dence of cosmic shear, the abundance of galaxy clusters, and lensing of the 
cosmic microwave background (CMB) photons. Observations of these quanti- 
ties can provide important constraints on growth, which in turn helps constrain 
cosmological models, and specifically distinguishes between dark energy and 
modified gravity. 

Notably, current constraints on the amplitude of mass fluctuations Sg and 
the function fog come largely from growth. They can be used to look for consis- 
tency of data with a cosmological model, as the Sg and fog(z) constraints are 
model-dependent. These constraints can also be used to look for internal con- 
sistency between different probes or datasets, as well as comparison between 
early- and late-universe constraints. Such comparisons have been very fruit- 
ful, and have led to the currently much-discussed Sg tension, where the CMB 
experiments appear to give a higher value of this quantity than cosmic-shear 
measurements. 

We highlighted several ways in which growth can be used to probe the 
cosmological model that go beyond simply measuring and reporting og, Sg or 
fos(z). Explicit parameterization of the growth sector is one such approach. 
For example, departures from growth predicted assuming general relativity 
can be enabled by freeing the growth index y from its general-relativistic value 
of ~ 0.55, or by parameterizing the gravitational potentials V and ® from 
their fiducial values. One can also encode additional freedom in the linear 
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growth function D(a) as a function of scale or time. Alternatively, one can 
explicitly separate information contained in the growth of structure from that 
in geometrical terms by introducing separate geometry and growth parameters 
in equations that govern cosmological observables. Because relations between 
the geometrical and growth quantities that are obeyed in LCDM and wCDM 
are likely to be broken in modified-gravity explanations of the accelerating 
universe, such “geometry-growth split” tests are potentially very powerful. 
Finally, we discussed how simply making predictions given by measurements 
by one set of probes (e.g. the CMB) for what another set of measurements 
(e.g. cosmic shear or some other, new probe, or an existing probe in a new 
range of scales or redshifts) is expected to measure, is also a very useful way 
to employ growth to stress-test the cosmological model. 

Most of the major upcoming experiments and telescopes with the role to 
probe the accelerating universe (see [156] for a recent review) rely on better 
measurements of growth, particularly at higher redshifts (z ~ 2-5) than cur- 
rently possible, in order to better understand dark energy and dark matter. 
Therefore, we expect that the measurements and interpretations of the growth 
of structure will become even more central to cosmology in the years to come. 
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